Direct observation of non-local effects in a superconductor 
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We have used the technique of low energy muon spin rotation to measure the local magnetic 
field profile B(z) beneath the surface of a lead film maintained in the Meissner state (z depth 
from the surface, z < 200 nm). The data unambiguously show that B(z) clearly deviates from an 
exponential law and represent the first direct, model independent proof for a non-local response in 
a superconductor. 
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Applying a small magnetic field parallel to the sur- 
face of a superconductor results in the expulsion of the 
magnetic flux from its interior, except for a small re- 
gion on the nm scale close to its surface where the local 
field B(z), measured at a depth z from the surface, is 
heavily damped (Meissner-Ochsenfeld effect). The spa- 
tial field dependence B(z) reflects the electromagnetic 
response of the superconductor and yields valuable infor- 
mation about its nature. If the Cooper pairs can be han- 
dled as point like, the electrodynamics of the system can 
be treated as local and the field extends exponentially 
over a typical length A (London penetration depth) of 
the order of a few tens to hundreds nm. Non-local effects 
in superconductors arise when the variation of the elec- 
tromagnetic field over the extent of the pairs cannot be 
neglected 0, . This is for instance the case in conven- 
tional superconductors with coherence length £ > A or at 
nodes of the energy gap of unconventional superconduc- 
tors, where the fc-dependent coherence length becomes 
effectively infinite [3j- In YBa2Cu3 06.95 indications of 
non-local/non- linear effects have been found in the field 
dependence of the effective magnetic penetration depth 
as determined byuSR measurements in the bulk of the 
vortex lattice 0, |5| • 

In conventional superconductors the penetration pro- 
file of the magnetic field in the Meissner state is predicted 
to clearly deviate from the usual exponential decay at the 
surface. Though theoretical predictions are already half a 
century old, the direct experimental verification has been 
lacking. Signatures for non-local field penetration have 
been searched so far by induction techniques [(| , magne- 
toabsorption resonance spectroscopy 0, Q , and polarized 
neutron scattering reflectrometry (PNR) 9]. The pio- 
neering measurements of Drangheid and Sommerhalder 
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[3 revealed that there is a sign reversal of B(z) (as pre- 
dicted by theory), but no further quantitative results 
could be drawn from this experiment. The magnetoab- 
sorption resonance spectroscopy technique uses the fact 
that quasi-particles travelling parallel to the shielding 
current are bound to the surface by an effective magnetic 
potential. Indication of non-local effects in Al was in- 
ferred by comparing microwave induced resonant transi- 
tions between the energy levels of these bound states with 
transition fields calculated from the energy levels of the 
trapping potential, parameterized to include the shape 
of the non-local BCS-like potential 0, @- Due to the 
resonant character of the experiment, only few specific 
points of the potential are probed. In addition the nor- 
mal metallic state has to be understood very well in order 
to interpret the data. This, together with uncertainties in 
modelling the surface bound states, makes the confirma- 
tion of the predicted functional form of B(z) rather indi- 
rect. The specular reflectivity of neutrons spin polarized 
parallel or anti-parallel to B also depends on the field 
profile. However, this technique requires model-fitting of 
spin-dependent scattering intensities rather than giving 
a direct measure of the spatial variation of the magnetic 
field. Up to now non-local corrections have been found 
to lie beyond the sensitivity of PNR Q . 

A direct measurement of B(z) requires experimental 
probes that allow to measure microscopically the mag- 
netic properties of the region extending only a few tens of 
nm away from the surface. We used the newly developed 
low energy muon spin rotation technique (LE-//SR) to 
map B(z) in superconducting lead [TfJ. We find that the 
functional dependence of B(z) is indeed non-exponential 
and that it follows the predicted Pippard-BCS theory 
HI S OH [IE [HI ■ We obtain a value for the London pen- 
etration depth Al = 57(2) nm and a clean-limit coherence 
length £ = 90(5) nm for Pb. 

A weak external magnetic field acts on the ground 
state of the superconductor as a perturbation. Within 
a perturbation expansion one can show [lil ITEl fl6l| that 
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the following non-local relation between the supercurrent 
density j and the vector potential A holds (in Coulomb 
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gauge V • A = 0): 
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where e is the electron charge, ns the supercarrier den- 
sity, to* the effective electron mass, and V A A = B. 
The first term in the square brackets, R a p, describes the 
paramagnetic response, whereas the second reflects the 
diamagnetic one. K a p is called the kernel. If the wave 
function of the electronic ground-state were "rigid" with 
respect to all perturbations (rather than only those which 
lead to transverse excitations) R a p would be identically 
zero and Eq.QJ would reduce to the local j-A relation 



(2) 



with fiQ the magnetic permeability of the vacuum. This 
combined with the Maxwell equation VAB = [igj yields, 
at a plane superconductor- vacuum interface, the result of 
an exponentially suppressed magnetic field 



B(z) = B cxt exp(-z/A L ) 



(3) 



with the London penetration depth Al = 

V Morris 

which is the well known result. 

However, R a p has a range of the order of the diame- 
ter of the Cooper pairs, i.e. of the coherence length £. 
The magnetic penetration depth sets the length scale for 
the decay of the magnetization; for A 3> £ the spatial 
variation of the vector potential A over the supercon- 
ducting pairs is negligible and the one-parameter local 
description of Eq.J2J holds. If £ > A the full non-local 
description has to be taken into account. Fourier analysis 
of the perturbation and of the response show that 
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B(z) = B c 



sin(gz) dq (4) 



where q is the wave vector, T the temperature, and £ 
the electron mean free path. The functional form of 
K(q£,T,£) (Fourier transform of the kernel in Eq.JU), 
starting from microscopic considerations, is explicitly 
known in the we ak ( BCS) 0, ^| and the strong cou- 
pling limit [Hill Il8| . The resulting formulae are rather 
involved but are veryclose to the phenomenological ex- 
pression of Pippard 0, [l5| 
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with x = q£(T, €) and 



J(0,T) 



1 



where £o 



€(T, t) £o V 
is the clean limit coherence length and J(0, T) is given 
according to Ref. 0. The kernel has the property 
fj,oK(x — > 0, T — * 0) = (corresponding to the local 
limit). This holds for the BCS kernel as well. 

Fig. Q]shows the theoretical predictions in the case that 
£ ^> A. The main features can be understood qualita- 
tively. In the non-local case the perturbing field changes 
over the extension of the Cooper pairs; since the charges 
within a Cooper pair do not experience the same force, 
the screening response is less effective and hence the field 
falls initially less rapidly than in the case of a point re- 
sponse. With the field penetrating further, at some range 
Cooper pairs "overcompensate" , which accounts for the 
curvature of \og(B(z)) as well as for the field reversal of 
B{z) Q. 

We have used the 100% polarized low energy muon 
beam at the Paul Scherrer Institute and the muon spin 
rotation technique (/xSR) to directly determine the val- 
ues of the magnetic field as a function of depth under- 
neath the surface [T^j. With a tunable energy between 
0.5 and 30 keV these particles are implanted one at a 
time at variable depth between ~ 1 nm and a few hun- 
dreds nm beneath the surface of the specimen. The local 
magnetic field B(z) at the stop position causes the muon 
spin to precess. The temporal evolution of the spin polar- 
ization of the muon ensemble, P(t), is monitored by the 
detection of the decay positrons which are anisotropically 
emitted preferentially in the direction of the muon spin 
at the moment of the decay 0] . This quantity is directly 
related by a Fourier transform to the internal magnetic 
field distribution sensed by the muon ensemble. The field 
distribution p(B) is connected to the implantation profile 
of the muons n(z, E) by: 



n(z,E) dz=p{B) dB 



(6) 



which states that the probability that a muon will ex- 
perience a field in the interval [B, B + dB] is given by 
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FIG. 1: Comparison between local and non-local magnetic 
penetration profiles for the same Al (parameters from Al for 
T -> 0, Al = 18 nm, £ = 1600 nm). The dashed-dotted line 
shows the typical exponential field profile predicted from the 
London theory. The solid line shows the non-local field profile 
from BCS theory. Specific features are: (i) The penetration 
profile is non exponential, (ii) The initial slope is less steep 
than in the local approximation, (iii) There is a field reversal 
before B(z) decays towards zero, (iv) The inset shows that 
log(_B(z)) vs. z, exhibits a clear curvature. 



the probability that it will stop at a depth in the range 
[z, z + dz]. Integrating Eq. JHJl on both sides yields: 
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FIG. 2: Measured magnetic penetration profile B(z) of Pb 
in the Meissner state for T = 3.05 K [T c = 7.21(2) K] in an 
external field of 8.82(6) mT, applied parallel to the surface 
of the film. Implantation energies E of the muons: 5.2 keV 
open circles, 14.8 keV closed circles. The mean values are 
plotted as closed squares with implantation energies (2.5, 4.0, 
5.4, 7.5, 10, 12.4, 15, 18, 22.5, 26, 30) keV from left to right. 
The solid curve is a fit according to Eq.l^J with the BCS 
parameters Al = 57(2) nm, £o = 90(5) nm, and I — 100 nm 
(fixed). The dashed curve shows the exponential dependence 
assuming the values from the full BCS fit. The inset displays 
results of Meissner state measurements of optimally doped 
YBa 2 Cu 3 7 -cS at T = 20 K (T c = 87.5 K). This curve is 
purely exponential, as expected. 
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which, for a chosen z, is an equation for B. We cal- 
culate n(z,E) with the Monte Carlo code TRIM.SP poj . 
which yields reliable implantation profiles for the muons, 
as shown in Ref. |2l|. Since n(z,E) is known and p(B) 
is measured, B(z) can be uniquely determined. 

To search for non-local effects we investigated thin 
films of type-I superconducting Pb. The films were sput- 
tered directly onto sapphire or quartz crystals mounted 
on a He flow cryostat. The samples had a diameter of 50 
mm and thickness of 1055(50) nm (sample I) and 430(20) 
nm (sample II) as determined by a high sensitivity sur- 
face profiler and Rutherford backscattering (RBS). For 
sample I an oxide-layer of 5.8(3) nm was found by RBS, 
whereas sample II had an oxide-layer of 16(2) nm. The 
critical temperature was determined by means of resistiv- 
ity and susceptibility measurements to T c = 7.21(1) K. 
The mean free path was estimated to be ~ 100 nm 
from resistivity measurements. After zero field cooling 
a -B ext = 8.82(6) mT was applied parallel to the sur- 
face. In the experiment low energy muons, with 
their spin perpendicular to the magnetic field and their 
momentum were implanted in the samples at variable 
depth up to 150 nm and the decay positrons detected 



by scintillation counters surrounding the sample. De- 
tails and characteristics of the low energy fi + source and 
spectrometer are given in p2|. The distribution p(B) 
was derived by maximum entropy Fourier analysis of the 
muon spinprecession frequency of the decay positron his- 
tograms m mm . 

The measured magnetic field penetration profiles B(z) 
in the Meissner state of Pb at T = 3.05(3) K and T 
= 6.66(3) K, respectively, are shown in the Fig. [3 and 
(sample I). Sample II gives consistent results. The 
curves show clear deviations from the exponential behav- 
ior with the characteristic curvature. Qualitatively, the 
initial slope of the curve is determined by A, whereas the 
log(B(z)) curvature is mainly governed by the ratio £/A. 
The deviations are more pronounced at lower temper- 
ature: on approaching T c the superconductor becomes 
more and more local. This is a consequence of the fact 
that A has a pronounced temperature dependence close 
to T c , but £ has not. The reversal of the penetrating field 
could not be detected since the muon range in the exper- 
iment was not large enough. We would like to point out 
that from Eq. Q it follows that the functional relation- 
ship B(z) can de determined from overlapping implan- 
tation profiles n(z, E) obtained at different energies. In 
order not to overload the figures only some energies are 
shown (displayed in different colors). Fig. 2 shows for in- 
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FIG. 3: Magnetic penetration profile B(z) of Pb in the Meiss- 
ner state for T = 6.66 K (T c = 7.21(2) K) in an external 
field of 8.82(6) mT, applied parallel to the surface of the film. 
Implantation energies of the muons: 14.8 keV open circles, 
30.0 keV closed circles. The mean values are shown as closed 
squares with implantation energies (2.5, 5.4, 7.5, 10, 15, 18, 
22.5, 26, 30) keV. The solid curve is a fit according to BCS 
theory, whereas the dashed curve shows an exponential with 
the same Al. 

stance curves obtained with 5.2 and 14.8 keV data. Their 
small relative deviation is a measure of the uncertainty in 
the derivation of B(z) and is an important sclfconsistency 
test of the reliability of the results. For a more detailed 
discussion and a compilation of additional data we refer 
to a forthcoming publication. Figs. |21 and also present 
the results from an alternative analysis using the mean 
values for (B) = J Bp{B) AB and (z) = J zn(z 7 E) dz. 
The points (B) plotted versus (z) are in very good agree- 
ment with the curve B(z) obtained from Eq.J7J). As a 



further cross check for the sensitivity and accuracy of 
the method based on Eq.Q we reanalyzed profiles in 
the Meissner state of optimally doped YBa2Cu30y_5 at 
T = 20 K (T c = 87.5 K, [^j). In this extreme type-II 
superconductor, at this temperature, the data perfectly 
follow an exponential law as expected for the local elec- 
trodynamic response (Fig|2| inset) and are also in very 
good agreement with the results previously obtained by 
an iterative solution of Eq.JBJ). Analyzing the Pb data 
with the BCS and Pippard theory leads to consistent re- 
sults for Al and £o ■ The strong electron-phonon coupling 
in Pb can be accounted for by a renormalization of the 
weak coupling parameters of the form Abcs - * / ^fZ 
and £bcs — ► £oZ with Z = 1 + A e - P h ~ 2.55 for Pb 
|27j (A e -ph, the electron phonon coupling). Compil- 
ing all the data, we find a magnetic penetration depth 
Al = 57(2) nm and a coherence length £o = 90(5) nm for 
T = K. 

In conclusion, by using spin polarized muons of a few 
keV energy as surface sensitive magnetic microprobes, 
we have shown that the magnetic penetration profile at 
the surface of superconducting Pb is non-exponential and 
that non-local electrodynamics effects, as predicted by 
Pippard and BCS theory, are responsible for this behav- 
ior. We believe that the ability to measure magnetic 
profiles beneath surfaces and buried interfaces on the nm 
scale, by means of LE-//SR, opens the door to explore 
interesting systems from fundamental as well as from ap- 
plied point of view. We look forward to further mea- 
surements of magnetic fields and fluctuations in surface 
superconductivity (e.g. proximity effects in multilayers, 
surface sheath) and magnetism. 
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[1] A. B. Pippard, Proc. Roy. Soc. (London) A216, 547 
(1953). 

[2] J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev. 

108, 1175 (1957). 
[3] I. Kosztin and A. J. Leggett, Phys. Rev. Lett. 79, 135 

(1997). 

[4] J. E. Sonier, J. H. Brewer, R. F. Kiefl, G. D. Morris, 
R. Miller, D. A. Bonn, J. Chakhalian, R. H. Heffner, 
W. N. Hardy, and R. Liang, Phys. Rev. Lett. 83, 4156 
(1999). 

[5] M. H. S. Amin, M. Franz, and I. Affleck, Phys. Rev. Lett. 

25, 5864 (2000). 
[6] K. E. Drangeid and R. Sommerhalder, Phys. Rev. Lett. 

8, 467 (1962). 

[7] R. E. Doezema, J. N. Huffaker, S. Whitmore, 
J. Slinkman, and W. E. Lawrence, Phys. Rev. Lett. 53, 
714 (1984). 



[8] R. E. Doezema, S. C. Whitmore, and J. N. Huffaker, 

Phys. Rev. B 34, 4614 (1986). 
[9] M. P. Nutley, A. T. Boothroyd, C. R. Staddon, D. M. 

Paul, and J. Penfold, Phys. Rev. B 49, 15789 (1994). 
[10] E. Morenzoni, F. Kottmann, D. Maden, B. Matthias, 

M. Meyberg, T. Prokscha, T. Wutzke, and U. Zimmer- 

mann, Phys. Rev. Lett. 72, 2793 (1994). 
[11] S. B. Nam, Phys. Rev. 156, 470 (1967). 
[12] S. B. Nam, Phys. Rev. 156, 487 (1967). 
[13] J. Halbritter, Z. Phys. 243, 201 (1971). 
[14] J. R. Schrieffer, Theory of Superconductivity (Addison 

Wesley, San Francisco, 1964). 
[15] M. Tinkham, Introduction to Superconductivity (Krieger, 

New York, 1980). 
[16] C. P. Poole, H. A. Farach, and R. J. Creswick, Supercon- 
ductivity (Academic Press, San Diego, 1995). 
[17] This expression is valid for specular reflection of the 



■5 



charge carriers at the surface, the case of diffuse scat- 
tering has been studied analytically and numerically and 
shows only minor differences for relevant parameters such 
as the penetration depth. 

[18] J. C. Swihart and W. Shaw, Physica 55, 678 (1971). 

[19] S. J. Blundell, Contemporary Physics 40, 175 (1999). 

[20] W. Eckstein, Computer Simulation of Ion-Solid Interac- 
tions (Springer, Berlin, 1991). 

[21] E. Morenzoni, H. Gliickler, T. Prokscha, R. Khasanov, 
H. Luetkens, M. Birke, E. M. Forgan, C. Niedermeyer, 
and M. Pleines, Nucl. Instrum. Methods 192, 254 (2002). 

[22] E. Morenzoni, H. Gliickler, T. Prokscha, H. P. Weber, 
E. M. Forgan, T. J. Jackson, H. Luetkens, C. Nieder- 
mayer, M. Pleines, M. Birke, et al., Physica B 289-290, 



653 (2000). 

[23] J. Skilling and R. K. Bryan, Mon. Not. R. astr. Soc. 211, 
111 (1984). 

[24] B. D. Rainford and G. J. Daniell, Hyperfine Interact. 87, 
1129 (1994). 

[25] T. Riseman and E. M. Forgan, Physica B 289-290, 718 
(2000). 

[26] T. J. Jackson, T. M. Riseman, E. M. Forgan, H. Gliickler, 
T. Prokscha, E. Morenzoni, M. Pleines, C. Niedermayer, 
G. Schatz, H. Luetkens, et al., Phys. Rev. Lett. 84, 4958 
(2000). 

[27] J. P. Carbotte, Rev. Mod. Phys. 62, 1027 (1990). 



